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QUALITATIVE  ASPECTS  OF  THE  MOTICSI  OF  RIGID  BODIES  WITH 
LIQUID-FILLED  TOROIDAL  CAVITIES 


ABSTRACT 


For  a  rigid  body  subject  to  no  momentB  the  differential  equations 


for  the  angular  velocity  can  be  solved  independently  of  the  remaining 
equations  of  motion.  The  integral  curves  are  Intersections  of  the  energy 


and  angular  m omentum  ellipsoids,  which  have  common  centers  and  principal 
axes.  In  general,  there  are  four  types  of  closed  integral  curves.  It 
is  well  known  that  if  the  solid  contains  a  cavity  that  is  topologically 
equivalent  to  (l.e.,  continuously  deformable  into)  the  interior  of  a 
sphere  and  Is  completely  filled  with  non-viscous  incompressible  fluid, 
the  properties  mentioned  above  remain  valid.  However,  If  the  cavity  is 
topologically  equivalent  to  the  interior  of  a  torus,  the  fact  that  the 
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new  possibilities.  The  angular  velocity  Integral  curves  are  still  inter¬ 
sections  of  ellipsoids  with  parallel  principal  axes,  but  one  of  the  centers 


has  been  displaced  through  a  distance  that  depends  on  the  parajneter  f*  . 

If  r  =  0  there  are  generally  four  types  of  closed  integral  curves;  five 
for  "weak"  circiilation;  three  for  "intermediate"  Ifl  ;  and  one  for  "strong" 
in.  The  qualitative  nature  of  the  Integral  curves  for  solids  with 


2; 


cavities  of  greater  topological  complexity  has  eilso  been  analyzed*  The 
number  of  distinct  types  of  behavior  is  surprisingly  limited  and  is,  in 
fact,  still  closely  akin  to  that  of  bodies  with  toroidal  cavities. 


k 


Page  intentionally  blank 


I 


Page  intentionally  blank 

h- 


Page  intentionally  blank 


1.  INTRODUCTION 


To  detenalne  the  motion  of  a  liquid-filled  solid  is  an  important  and 
difficult  problem  in  dynamics.  In  actual  applications  the  liquid  is 
viscous  and  need  not  completely  fill  the  cavity  in  which  it  is  stored. 
However,  when  one  attempts  to  attack  this  problem,  a  natural  model  with 
which  to  begin  is  that  of  a  cavity  completely  filled  with  an  inccanpressible, 
non-viscous  liqxxld  in  irrotatlonal  motion. 

If  the  cavity  is  topologically  spherical,  (i.e.,  homeomorphic  to  the 
Interior  of  a  sphere),  then  in  this  model  the  liquid  filling  introduces  no 
dynamic  novelties.  The  differential  equations  of  motion  of  the  solid- liquid 
system  can  simply  be  considered  to  be  those  of  a  solid  with  a  different 
mass  and  different  inertial  matrix.  Among  its  other  shortcomings  there  is 
no  possibility  for  including  in  this  unsteady  irrotatlonal  flow  a  "spin" 
or  general  circulation,  such  as  one  would  expect  to  obtain  if,  for  example, 
the  container  itself  is  spinning. 

If  the  cavity  is  topologically  equivalent  to  the  interior  of  a  torus, 
it  becomes  possible  to  have  unsteady  irrotatlonal  flow  with  circulation, 
r  .  Unfortunately  T  must  be  constant,  a  fact  which  seriously  restricts 
prospective  application  of  the  following  theory.  Nevertheless  it  seems 
worthwhile  to  Investigate  the  influence  of  circulation  on  the  motion  of  the 
composite  solid-liquid  system.  For  this  purpose,  the  most  natural  and 
easiest  problem  to  consider  is  the  classical  one  of  motion  under  no  forces 
and  moments . 

If  we  consider  a  system  with  no  special  symmetry,  the  results  may  be 

2 

summarized  as  follows.  As  will  be  shown  later,  the  range  of  T  must  be 

2  2 

divided  into  three  intervals  with  end  points  0,  *~2  ’  *~1  ’ 

r  =  0  we  obtain  the  classical  result  that  there  are  six  possible  steady 

states  of  rotation,  four  of  them  stable,  the  other  two  mstable.  Roughly 

speaking,  there  are  four  periodic  types  of  angular  motion,  each  centered 

about  one  of  the  stable  steady  states.  For  weak  circulation, 

2  2 

0  <  r  <  ^2  *  "tbere  are  still  six  steady  states,  four  of  them  stable. 
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Now,  however,  in  addition  to  four  previously  mentioned  types  of  periodic 

motion,  each  centered  about  one  of  the  stable  steady  states,  there  is  now 

a  fifth  type  of  periodic  motion  centered  about  two  stable  steady  states. 

2  2  2 

For  intermediate  circulation,  ^^1^  there  are  four  steady 

states,  three  of  them  stable,  and  three  types  of  periodic  motion.  Finally, 

2  2 

for  strong  circulation,  <  F  ,  there  are  two  stable  steady  states, 

and  one  type  of  periodic  motion. 

Certain  types  of  degeneracy  are  associated  with  appearances  of 
multiple  characteristic  roots  of  a  matrix  involved  in  the  theory,  or 
special  orientations  of  an  "axis"  vector  associated  with  the  cavity.  An 
important  example  of  such  degeneracy  occurs  for  axlsymmetrlc  liquid-solid 
systems . 

If  the  cavity  Is  of  topological  genus  N,  i.e.  topologically  equivalent 
to  the  interior  of  a  sphere  with  N  handles,  the  analysis  of  the  possible 
motions  can  be  performed  by  an  extension  of  the  methods  developed  for  the 
discussion  of  toroidal  cavities  (N  =1).  With  a  cavity  of  genus  N  we 
associate  N  "axis"  vectors,  of  which  not  more  than  three  can  be  linearly 
Independent,  of  course.  This  yields  the  conclusion  that  there  are 
essentially  four  types  of  motion  for  bodies  with  liquid-filled  cavities  of 
genus  N,  depending  on  the  number,  between  zero  and  three,  of  independent 
"axis"  vectors.  Two  types  correspond  to  topologically  spherical  or  toroidal 
cavities  (N  =  0  or  l) .  We  conjecture,  without  proof,  that  the  other  two 
types  correspond  to  cavities  of  genera  N  =  5  or  4.  The  measure  of  arbi¬ 
trariness  available  for  the  choice  of  circulations  for  cavities  of  genus 
N  =  4  merely  serves  to  provide  a  greater  variety  of  assignments  of 

circulation  that  will  produce  the  same  motions  as  a  body  with  a  cavity  of 
genus  N  ^  3. 
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2.  FLOWS  WITH  MOVING  BOUNDARIES 


The  most  commonly  discussed  problem  that  involves  rigid  bodies  and 

a  non-vlscous  incompressible  fluid  concerns  the  motion  of  a  finite  solid 

(2  5’) 

throxigh  an  infinite  region  filled  with  fluid'  ’  ' .  The  method  used  for 
this  classical  problem  requires  only  slight  modification  to  adapt  it  to 
develop  the  equations  of  motion  of  a  finite  rigid  body  that  contains  a 
cavity  completely  filled  with  fluid.  However,  since  our  pvirpose  in  the 
following  discussion  is  to  emphasize  the  influence  of  circulation  and 
since  the  derivation  is  comparatively  short,  for  the  sake  of  con5)leteness 
and  clarity  we  shall  reproduce  it  here. 

First  choose  a  moving  system  of  rectangiilar  coordinates  relative  to 
axes  rigidly  attached  to  the  body  and  with  origin  at  the  center  of  mass 
of  the  combined  mass  of  body  and  fluid.  At  any  time  t  also  choose  coordi¬ 
nates  (x,y,z)  relative  to  a  system  of  axes  fixed  in  space  that  coincides 
with  the  Instantaneous  position  of  the  moving  axes.  At  this  instant  the 
position  of  any  point  in  the  body  or  the  fluid  can  be  specified  by  its 

coordinate  vector  x  =  (x,y,z).  Also  at  time  t  let  u(t)  =  (u  ,  u  ,  u  ) 

~  X  '  y  z 

be  the  velocity  of  the  moving  origin,  and  let  w(t)  =  (w  ,  w  ,  w  )  be  the 

*—  X  y  z 

angular  velocity  of  the  moving  system  relative  to  the  fixed  axes. 

Wow  suppose  the  body  contains  a  cavity  filled  with  incompressible 
nonviscous  fluid  of  density  Let  be  volume  occupied  by  the  cavity, 
and  its  boundary.  If  we  assume  that  the  motion  of  the  liquid  was 
sta2rted  impulsively  from  rest,  then  by  Kelvin's  theorem  it  will  be 
irrotational.  Accordingly,  the  velocity  of  the  liquid  can  be  eDopressed 
as 


L  — ^  ~  V  -ar  \,x,  u; 

for  some  velocity-potential  function  such  that 

V  =  0  '  in 

Ij 


(2.1) 


(2.2) 
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and  subject  to  the  condition  that  on.  the  normal  component  of  the 
velocity  of  the  fluid  relative  to  that  of  the  boundary  vanishes.  Let 
X_  be  any  point  of  and  n  the  corresponding  unit  normal  directed  into 
V_ .  Then  we  must  have 

XJ 


u^/bn  =  n  =  V  v|>  (x,  t)  =  n  * 


bince  n  =  n(x)  at  time  t,  tentatively  choose 


u^^t)  +  v(t)  XX  on  S3  .  (2.3) 


,  t)  =  u( 


J.\  .  J  f—\ 

b;  V"  • 


x)  4-  w(t)  ‘  cr(x) 


Then  (2.2)  and  (2.3)  will  be  satisfied  if  we  choose 
single-valued  solutions  of 

^  =  0,  b  =  0  in  V. 


and 


3  <i>/bn  =  n  •  V  ^  =  n 
5  =  n’Vu  =  xxn 


L 


on  S, 


(2.t) 

(x)  and  a  (x)  to  be 

(2.5) 

(2.6) 


If  ^  is  topologically  equivalent  to  the  interior  of  a  sphere,  then 
<!>)  and  cr  must  all  be  single -valued,  and  as  solutions  of  Meumann  problems 
they  must  be  unique  except  for  additive  constants.  If,  however,  V  is 
topologically  equivalent  to  the  interior  of  a  torus,  and  if  we  overlook 
the  question  how  one  would  create  a  general  circulation  in  the  cavity, 
and  thus  abandon  the  impulsive  start  from  rest,  then  need  no  longer 
be  single -valued.  This  can  be  seen  by  considering  the  circulation 


r~  r 

' 


C)  =  j  (x,  t) 


ox  = 


r 

J 


a  ^ 


{2.f) 


wVif 


1^“  «'j  "n  T  r 

XU 


o  .T  1  r  I  C 
oxiup-nj  — (ji_UD 


cu  path  C  li-i  v'  .  Now  form  r  (C‘)  for  any  other 

Xj 


simply  closed  path  C '  which  can  be  continuously  deformed  into  C  without 
crossing  Sj^.  If  'under  this  deformation  the  sense  in  which  C '  is  traversed 


in  r  (n 


spends 


sense  of  C  in  r  (C),  then  by  Stokes'  theorem 


r  (u')  =  r(C);  otherwise  r  (C)  =  -  r  (c).  If,  in  particular  C  can 
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be  continuously  deformed  into  a  point  within  V  ,  then  r  (C)  =  0.  On 
the  other  hand,  if  C  loops  once,  and  C  loops  N  times  in  the  same  sense 
about  the  hole  of  the  torus,  then  r  (C')  =  N  r  (C).  Thus,  for  all 
single  loops  traversed  in  the  same  sense  we  get  the  same  circulation 
r(C)  =  r  .  If  r  ^o,  (2.7)  implies  that  <I>  must  be  multiple -valued. 
Accordingly,  let  us  modify  (2.4)  to  the  form 

(x,  t)  =  u(t)  •  ^(x)  +  w(t)  .  ff(x)  +  r  t(x)  (2.8) 

where  in  addition  to  (2.5)  and  (2.6)  we  require 

V  =  0  in  V  (2.9) 

Xj 

br/ba  =  n  •  V  t  =  0  (2*10) 

and  though  the  ccanponents  of  V  t  are  single-valued,  t  increases  by  imity 
when  a  closed  loop  about  the  hole  of  the  torus  is  traversed  once  in  an 
arbitrarily  selected  positive  sense.  By  Bernoulli's  theorem 

-  p/pj^  =  8<I)/8t  +  0.5  (V<I>)^  +  ^(t)  •  X  +  F(t)  (2.11) 

where  G(t)  is  a  unit  vector  parallel  to  the  gravitational  field  and  P(t)  is 
some  function  of  t.  Since  the  pressure,  p,  must  be  single-valued  in  7,, 

Xj 

then  d  r  /dt  =  0,  l.e.  T  must  be  a  constant . 

Hereafter  it  will  suffice  to  consider  (2.8)  without  special  reference 
to  (2.4),  since  results  for  topologically  spherical  cavities  can  be  deduced 
by  merely  setting  F  =  0  in  the  following  discussion. 

For  general  cavities  Vj^  of  finite  extent,  which  are  not  necessarily 
even  topologically  toroidal, 

/  (x)  =  X  (2.12) 

satisfies  (2.5)  and  (2.6).  The  discussion  of  £(x)  and  t(x),  however, 
cannot  be  continued  without  specializing  Vt*  will  be  done  for 

jj 

axisymmetric  cavities  in  Section  7* 
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5.  LINEAR  AND  ANGULAR  MO^TA  OF  THE  LIQUID 

To  formulate  the  equations  of  motion  of  our  composite  solid-liquid 
system  we  shall  require  the  linear  and  angular  mcsmeata  of  the  liquid. 

Let  the  linear  momentiam,  and  k  an  arbitrary  constant  vector.  Then 


11  ■  -  -  f  V  ^  +  J2  '^3 


(3.1) 


where ,  in  accordance  with  (2.8)  and  (2.12),  the  scalars  J  are  defined 
below.  First 


P^-k  •  u  dV  =  M,  k  •  u 
’  ’L  ^  ~  ^  - 


¥  » 


Thus 


Jp  =  Mr  *  2St  ° 


where  3^  is  the  center  of  mass  of  the  liquid.  Finally 


f5.2^ 


where  is  the  total  mass  of  the  liquid.  Next 

'^2'^^L  ""  /v  -  *  ^  (w  •  c)  dV  V  (k  •  x)  .  V  (w  •  a)  dV 

L  L 

By  Green's  theorem  and  (2.6) 

'^2/^L  “  "  f  S  -  2)/^ 

Xi 

=  -  /  „  (k  •  x)  (w  X  X  •  n)  dS 
j  Sl-  -  -  -  - 

=  J  y  ”  ’  [  (±  '  i  J 


Jjp. 

J  n 


fv  r-k.V  Tdv  =  rf  V  (k  •  x)  •  7  V  dv 

»  ^T.  ¥  ’t  “■ 


=  -rf 


(k  •  x)  dS 


10 


whence  by  (2.10) 


J.  =  0 

D 


(5A) 


Since  k  was  arbitrary,  these  results  imply 


t 


T  W  JL 


\J-^I 


Next,  let  Tj-(t)  be  the  angiilar  momentum  of  the  liquid.  Then 

Ti  =  ^  D  X  X  dy  =  -T  +  .T  +  .T 

-^L  J  ^  ■  -5  ■  ^ 

li 


for  the  following  choices  of  vectors  J.  First 

h  =  /vPl^^ii‘iV  =  Mi^xu  (5.7) 

Jj 

Next,  by  Green's  theorem  and  (2 *6) 

J^/Pl  =  Jy  X  X  V  (w  •  ff)  dV  =  -  J*  y  V  X  1^  (w  *  a)  xl 

L  Xj 

=  (w  •  a)  X  X  n  dS  =  -  (w  •  o)  dS 

(5.8) 


dV 


=  V  (w  '  ff)  •  V  d  <iy 

J 


Finally. 


f  =  /„xxVT'<iV  =  -  /^Vx  (tx)  dV 

^  X.  ^  -L  ^ 

By  means  of  some  surface  3*  bounded  by  a  closed  curve  C*  on  3^  change  V 

1j  l 

into  a  topologically  spherical  region  in  which  t  is  single  valued.  Let 
n^  and  denote  the  inward  unit  normal  and  value  of  t  on  the  "initial" 
side  of  S*',  and  let  +  1  be  the  corresponding  functions 

at  the  same  points  on  the  opposite  or  "final" side  of  3^.  Then 

Jg/Pl  r  «  J*/Pl  =  ~  fs  (x  X  n)  T  dS  +  J*  g^(x  X  n^)  dS 

’  ’  (3.9) 
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Later,  in  order  to  simplify  the  equations  of  motion,  we  shall  make 

transformations  of  coordinates  that  can  be  most  easily  motivated  by 

exploiting  the  relation  between  ^  ,  n  ,  and  the  kinetic  energy  T  of 
(2  — L  -^L  L 

the  liquid'  ’  ' .  We  have 


1j 


(3.10) 


where 


=  f  Y 


dV 


(3.11) 


Jo  =  2  r  ru  +  ^  w  .  ct)  •  V  T  dV 


y  Vl 


By  Green ' s  theorem 


_  2  /*  ,  n_\2  _  2  _  * 

(  J  V,  ^  ^  ‘^9 


=  2p  r  r  „  V  (u  •  X  +  w  •  ff)  .  V  T  c 

^  ^  'L  -  -  - 

=  -  2p^  r  (u  .  X  +  w  •  ff)  Sx/Bn 


whence  by  (2.10) 


dS 


(3.12; 


T.  =  n 


f  Ti  -I  ■Z\ 


Now  observe  that 


St, /3u  = 

h'  — 


r_.  o-  y  <J)  dv  =  e 

J  V^L  •  -  • 


Also 


1  )i  ^ 

V, 


STt./Sw  =  f  V4>  •  V  a  dv  =  r  D_  (u  +y  w  •  fT  +  rUx;  .  y  «  (iv 
n-  J  Y  Ij  -  -  J  YjL  ^ 


If  in  Jg  we  set  u  =  0  we  obtain  for  any  w 


J  ^  Pj.  rVT-Vff  dV  =  0 

-  -L 
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Thus  by  Green's  theorem  and  (2.6) 


Hence 


(dTj^/2rw)/pj^ 


=  y*„V(u*x  +  w*£)*V£  dv 

'  ’L 

=5  -  r  „  (u  •  X  +  w  •  c)  X  X  n  dS 

;Sl-  -  - - - 

=  -  r  „  V  X  I  (^  *  X  +  w  •  a)  X  I  dV 

j  L  -  -  J 

=  /*  X  X  V  (  <I>  -  r  t)  dV 

'  'L 


=  ar/atw  +  r  £»  (5.15) 

where  we  shall  call  J*  the  axis -vector  associated  with  t. 

It  will  prove  convenient  to  write  T^^  in  matrix  notation.  First 
observe  that 


^  +  2  J 


v/l  ^ 

jj 


7  w 


dV  + 


0 

J 


(V  w 


of 


dV 


Thus,  if  we  interpret  u  and  w  as  column  vectors,  then 

2Tj.  =  Mju"^  lu  +  2^^  w  +  w'^Cj^  w  +  J5*  (5.l6) 

where  I  is  the  5x3  identity  matrix,  and  are  matrices  with  constant 
elements ,  and 


where  the  superscript  T  denotes  "transpose",  and  J^*  is  a  constant.  Since, 
by  repeatedly  used  manipulations 


15 


u^Bt  w/p^  =  /*u  ■  V  w  •  a  dV  =  -  f  (u  ■  x)  b  (y  •  a)/bn  dS 


/s  ^  £ 


dS 


•  f  V  X  (u  •  x) 


-  J  \  ' 


I  V  (£  ■  H  X  x)  dV 

^  'L 


X  dV 


T>  TJ 

"L- 


M  U  •  1 '  "v  V 

‘•*1  £  h  -  ^ 


f  ^  1 

J-wy 


Ilf 


k.  LINEAR  AND  ANGULAR  MOMENTA  OF  THE  SOLID 


Let  be  the  volijme  occupied  by  a  solid  of  mass  Mg,  density  pg, 

and  with  center  of  mass  at  x„.  Then  the  linear  mcanentum  of  the  solid  is 

—a 


-S  ""  fv  Pg  (ii  +  V  *  *)  ‘iV  =  Mg  (u  +  w  X  2^) 


(V.l) 


and  its  angular  momentum  is 


T>_  =  r  o_  X  X  (u  +  w  X  x)  dV 


C  r  2  ,  V  1  ... 

U  +  J  V  I  -  -  "  -  J 


(4.2; 


Its  kinetic  energy  T^,  cein  be  calculated  from 


2T. 


=  f  y  pg  (u  +  w  X  x)^  dV 


=  MgU^  +  2MgU  •  w  X  3^  +  J  „  pg  -  (x  .  w)‘'l  dV  (4.5) 

s 


Clearly 


=  ar^/au 

-nD  O  — 


=  dTg/aw 


(4.4) 

(4.5) 


In  matrix  notation  T  takes  the  form 


2T^  =  Mq  u'^  Iu  +  2u'^B_  w  +  w"C  W 

li—  —  — •S— •  ~S— • 


(4.6) 


where  and  C  are  constant  matrices,  and 

s  s 

p  T  _  ^ 

“S  ’S 

T 

u  B_  w  =  M_  u  •  w  X  x_ 

—  S  —  D  —  —  — S 


(4.7) 

(4.8) 
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5.  EQUATIONS  OF  MOTION 


For  the  composite  liqiild- solid  system  the  total  kinetic  energy 


T  =  0.5  Mu^  lu 


+ 


T 

w  +  0.5  w  c  w  +  r 


(5.1) 


where  M  =  M,.  +  M  is  the  total  mass,  6=3^+  B_,  and  0=0,.  +  C_  are 

1j  o  1^0  Li  o 

constant  vectors,  and 


_T 

0 


0 


(5.2) 


Since  T  must  be  a  positive  definite  quadratic  form  in  r  and  the  components 
of  u  and  w,  then  as  a  matter  of  fact  0  must  also  be  positive  definite 
(3.18)  and  (4-. 8)  imply 


T_ 

u  nw 


“ 

u  •  W  X 


''sis. 


(5.3) 


Hereafter  we  shall  assume  that  the  origin  is  at  the  center  of  mass  of  the 
composite  system.  Thus  M^x..  +  M^.,  =  0,  which  implies 

U—Lt  CJ—O 


B  =  0  (5.4) 

The  linear  and  angular  momenta  are  1  =  ii  “  jlL  iJs’  ^ 

accordance  with  the  results  of  the  preceding  sections 


i  =  ar/Bu  =  Mu  (5.5) 

T1  =  ar/aw  +  r  J*  =  Cw  +  rj*  (5.6) 


If  F  and  L  are  the  resultants  of  the  external  forces  and  moments  acting  on 
the  composite  SyS'fcSiu  j  then  the  equations  of  motion  become 

M(u’  +  w  X  u)  =  F  (5.7) 

c  w'  +  w  X  (c  w  +  r  j*)  =  L  (5.8) 

where  '  denotes  d/dt.  Note  that  C  is  entirely  determined  by  the  geometry  of 
the  system  and  its  mass  distribution.  The  circulation  of  the  fluid  manifests 
itself  only  in  the  term  F  J*  of  (5.8). 


6.  MOTION  SUBJECT  TO  NO  EXTERNAL  TORQUE 


The  motion  of  a  rigid  body  in  the  absence  of  external  torques  is  a 
standard  topic  for  mechanics  textbooks ' .  Accordingly,  the  case  L  =  0, 
which  oeeur-s  for  example  for  motion  in  a  uniform  gravitational  field, 
should  be  ideally  suited  to  bring  out  very  clearly  the  novelties  introduced 
by  the  Inclusion  of  circulation. 

For  L  =  0  (5.8)  becomes 


CW+wx(Cw  +  rJ*)  =  0 


(6.1) 


Accordingly  it  is  possible  to  determine  the  angular  velocity  w(t)  indepen¬ 
dently  of  u(t).  This  system  has  two  well-known  integrals 

w"C  w  =  2‘J^  (6.2) 

and 

(c  w  +  r  j*)'^  (c  w  +  r  =  K  (6.3) 


where  the  constant  T*  is  that  part  of  the  kinetic  energy  associated  with 
w,  and  the  constant  K  is  the  square  of  the  magnitude  of  the  angular 
momentiim  vector.  If  the  rectangular  components  of  w  are  interpreted  as 
rectangular  coordinates,  (6.2)  is  an  ellipsoid  with  center  at  w  =  0,  and 
(6.3)  is  an  ellipsoid  with  center  at  w  =  -  T  C’'^J*  which  is  in  general 
different  from  w  =  0.  Thus  the  integral  curves  of  (6.1)  are  intersections 
of  ellipsoids,  but  by  contrast  with  the  more  familiar  case  r  =  0,  the 


B  of  the  ellipsoids  no  longer  coincide. 


It  should  be  remarked  that  under  the  transformation  w(t)  =rfi(rt) 


utiA-co  out?  iOna 


Cdfi/a  rt  +  X  (Cfi  +  J*)  =0 

Thus,  if  the  structure  of  the  solutions  of  (6.1)  has  been  determined  for 
one  value  of  •“  /  0  ,  then  it  is  known  for  all  r  ^  0.  We  shall  not  exploit 
this  fact,  however.  In  the  sequel. 
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Hereafter,  let  us  suppose  that  has  been  prescribed.  To  gain  a 
comprehensive  viev  of  the  behavior  of  the  associated  one-parameter  family 
of  integral  curves  for  given  r ,  let  us  consider  the  level  curves  of  the 
function  K(w),  defined  by  the  left  member  of  (6.3),  on  the  energy  ellipsoid 
(6.2).  Let  us  begin  by  determining  the  critical  points  on  (6.2)  for  which 
K(w)  is  stationary.  Proceed  by  Lagrange's  method  of  undetermined  multipliers. 
Let 

f(w,  x)  =  (c  w  +  rcp)"^  (c  w  +  r  j*)  -  x(w'^cw  -  21^) 


Then  at  the  desired  critical  points 

ap/5iw  =  20*^  (Cw  +  r  J*)  -  2X  C  w  =  0 

rp 

and  Sf/Sx  =  0,  which  merely  reasserts  (6.2).  Now  by  (5.2)  C  =  C.  Since, 
furthermore,  C  is  non-singular,  then  SF/dw  =  0  Implies 

(C  -  X  I)  w  =  -  r  J*  (6.4) 

If  X  is  not  a  characteristic  root  of  C,  then  the  critical  points  of  K(w) 
are  defined  by 

^  =  -  (C  -  X  1)“^  r  J*  (6.5) 

as  functions  of  X.  Of  course  these  are  merely  the  singular  points  of  the 
system  (6.1).  Then  (6.2)  and  (6.5)  imply 

(c  -  X  l)"^  c(c  -  X  1)“^  J*  =  2T*/r  ^  (6.6) 

2 

from  which  X  must  be  determined.  In  accordance  with  (5-1)  T^/  r  is 
proportional  to  the  ratio  of  the  kinetic  energy  due  to  w  to  that  due  to  r 

To  simplify  our  discussion,  let  us  assume  that  the  x'  coordinate 
axes  have  been  chosen  to  be  parallel  to  the  principal  axes  of  C.  Then 


C  = 


and,  in  general,  we  may  assime 


0  <  D  <  E  <  F 


(6.7) 
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T 

If  we  also  assume  J*  =  '(j>  k,  i),  then  (6.5)  and  (6.6)  take  the  special 
forms 

w  =  (j/x-D,  k/x-E,  i/X“P)  (6.8) 


and 

Dj^(X-D)"^  +  Ek^  (X-E)“^  +  f/  (X-F)“^  =  21^/ T  ^  (6.9) 

2 

If  in  (6.9)  ve  consider  T*/  r  to  be  a  function  of  X,  then  since 
2  ^2  2  -*2 

d  (G?^  r  )/dX  >0  always,  d(T^  T  "  )/dX  always  increases.  Accordingly 

_2 

the  graph  of  I*  f”  has  the  qualitative  form  shown  in  Figure  1,  with  the 

_2 

horizontal  asymptote  1*  T  =0  and  three  vertical  asymptotes  X  =  D,  E, 

-2  -2  -2 
or  F.  Let  T  and  T*'  r  be  relative  minima  of  F  ,  and  assume 

-2  "^-2  ^ 

r  -  <  1*  F  ^  ,  thotigh  the  sense  of  this  ineqmlity  is  not  essential 

^  c.  2 

in  the  sequel.  Then,  depending  on  the  size  of  F  ,  (6.9)  will  have 

from  two  to  six  real  roots.  From  Figure  1  it  is  apparent  that  the  smallest 

real  root  will  always  be  less  than  D,  and  the  largest  always  exceeds  F. 

2 

In  the  limiting  case  F  =  0,  we  obtain  three  double  roots  X  =  D,  E,  and 
F,  which  by  (6.4)  correspond  to  the  ends  of  the  principal  axes  of  the 
energy  ellipsoid  (6.2),  as  indicated  in  Figure  1. 


The  general  nature  of  the  stationary  values  K(w  )  caii  be  determined 


as  follows.  Let  N  ,  N  ,  and  N  ,  respectively  be  the  number  of  relative 
M  s' 

maxima,  relative  minima,  and  saddle  points  of  K  on  (6.2).  If  these 

(4) 

numbers  are  finite,  the  Morse  theory  of  critical  points  asserts 


N 


m 


^  1 


N, 


M 


^  1 


N  -  N  =1, 
ms 


Nw  -  N 

M  s 


(6.10) 


m  M 


N 


=  2 


If  we  disregard  occurrences  of  double  roots  X,  which  will  be  taken  into 
account  later,  (6.10)  yields  the  possibilities  tabulated  hereafter: 
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Criulcul  ^joinus 


M  in 
N 


u 


k 

2 


II 

I. 

5 

1 


III 

2 

0 


When  there  are  only  two  critical  points,  of  course  N  =  N.,  =  1. 

m  M 

In  the  classical  cases  of  rigid  body  motion  or  of  motion  of  liquid 

xj.-Lxc:;u.  u^^u.j.c^c;>  wxuiiv4u.u  uxx  i^ujucx.  iv  licto  owl/  lucLjLXilia.  cxu  one;  c:;iiu.&  l/±  uiicr 

minor  axis  of  the  energy  ellipsoid,  two  minima  at  the  ends  of  the  major 
axis,  and  two  saddle-points  at  the  ends  of  the  mean  axis.  K  has  the 


ai*  ‘Krv'h'h  gq<^<^1<=»  T^r\•^r^4‘G  oyi<^  +Vn: 

w  ky  W  OJ.A  -  j.y J.  ^  ‘ 


r'livnr^i  — 

iV  — 


the  surface  of  the  energy  ellipsoid  into  four  parts,  each  of  which 
contains  a  family  of  closed  Integral  (or  level)  curves  surrounding  one 
of  the  maxima  or  minima.  Figure  2a  is  a  schematic  representation  of  the 
system  of  integral  curves  on  a  cut  and  flattened  ellipsoid. 


As  r  varies  continuously  the  locations  of  the  critical  points  will 

Continuous  depen- 


<22 

vary  continuously  on  (6.2)  as  long  as  0  =1”  ^ 


dence  on  F  will  assure  that  maxima  of  K  move  into  maxima,  minima  into 

2 

minima,  and  saddle-points  into  s addle -points .  For  small  F  >  0,  however 

^  ^  T  ’  1  ^  V  4”  4*  Vi  /“V  4*t  n  ^  ■»*\y^4  4*  4-  4  01/—.  ynT  »  4»  V.  4  n 

ullc;  vcva.U.co  v/i.  i\  cv  u  oiic:  owv/  oa.vj.<uUi.c  mUo  u  •  j.w  biiww 

observe  that  by  (6.3)  and  (6.4)  we  have 


o  o 

K  =  X"w" 


=  K* 


at  critical  points,  and  then  by  (6.8)  and  (6.9) 

.■  a-£)  -A.  -  r^-i^ 


(X-D)"" 


(X-F) 


^  J 


Then 


r  Omit 


^  =  2x1-^  .(l_|)_JJi^  +(|.1)_A^| 

^  (X-D)^  ^  (X-F)^J 


(6.11) 


(X-D)' 


(4-F)' 
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_  2  _  2  _  2. 

Let  2e  =  min  (E-D,  F-E).  Then  for  some  Fe  min  (  >  ^2  ^ 

shall  have  X  "^dK^/dX  >  0  for  all  T  r  e  and  |  X-E  |  ^  6.  Furthermore, 

2  <  2 

there  exists  some  r  =  Te  such  that  (6.9)  will  have  exactly  two 

o  2  2 

roots  in  |X-E|<£  for  r  <  T  .  Since  dK^/dX  >  0  on  the  interval 

o  2  2 

joining  these  roots,  this  implies  that  for  0  <  r  <  r  ^  the  values 
of  K*,  say  and  are  different. 

Now  the  level  curves  K  =  K,  and  K  =  /  K,  must  continue  to 

X  d  X 

resemble  lemnlscates,  with  double  points  at  the  saddle  points.  Since 

_  2 

they  cannot  intersect,  the  situation  for  small  r  must  resemble  that 
shown  in  Figure  2b.  Th^  four  families  of  integral  curves  surrounding  a 
maximum  or  minimum  within  one  of  the  lobes  of  K  =  or  are  obviously 
counterparts  of  families  encountered  for  f“  =  0.  'Pne  novelty  introduced 
for  small  r  ^  >  0  is  the  occvirrence  of  a  fifth  set  of  closed  Integral 
curves,  typified  by  the  dashed  cixrve  between  the  level  curves  K  = 
and  K  =  K^.  Although  our  proof  that  is  valid  only  for  sufficiently 

small  values  of  T  it  seems  plausible  that  the  result  is  true  for 

2  2 

0  <  r  <  r  2  . 

2  2 

As  suggested  by  Figxire  1,  when  T  ~  ^2  maxima  (for 

the  conditions  depicted  in  our  graph)  should  coalesce  with  the  saddle- 

'2  2  2 

point  Sg*  For  f”  g  ^  ^  *~1  ’  'there  remain  two  minima,  one  maximum, 

HTK^  one  saddle-point.  Now  there  will  be  three  types  of  closed  integral 
2  2 

curves.  When  T  =  P  -■  one  of  the  minima  will  coalesce  with  the 

i  2  2 

remaining  saddle-point.  Finally,  for  P  ^  <  T  ,  there  remain  one 
maximum  and  one  minimum.  Now  there  is  only  one  set  of  closed  integral 
curves.  From  Flgvure  1  it  is  clear  that  |  X|  and  |r|  tend  to  infinity 
together .  Since  for  large  |  X |  we  can  expand 


^  n 


V'-' 


1 


/V 


then  by  (6.6) 


p  p 

2T^<-X  /  P 
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llm  {  r /\)‘^  =  2T*-/j*^CJ* 

f r(=  oo 

Then  by  (6.5) 

lim  V  =  +  (2T^/j*'^CJ^)°‘^J* 

iri=  oo 

i.e.,  bhe  critical  points  for  j  Ti  =  cx>  are  at  the  ends  of  the  diameter  of 
the  energy  ellipsoid  parallel  to 

In  the  classical  case,  T  =  0,  the  components  of  v  can  be  expressed 
in  terms  of  elliptic  functions.  For  r  /  0  and  D,  E,  F  all  distinct  this 
no  longer  appears  to  be  the  case.  For,  let  us  attempt  to  express  w  and 
w  as  functions  of  w  on  an  integral  curve.  When  we  eliminate  v  ,  for 

Z  X  z 

example  from  (6.2)  and  (6,3),  we  shall  obtain,  in  general  an  equation 
involving  a  polynomial  of  fourth  degree  in  w  ,  in  which  all  powers  of 

y 

w  betveen  zero  and  four  can  actually  occur •  Thus  v  becones  a  complicated 

y  y 

algebraic  function  of  w  ^  and  so^  presumably,  does  v  .  Thus  integration 

X  z 

of  the  equation  of  motion  that  involves  w  will  presumably  lead  to  some- 

^  V 

thing  more  complicated  than  elliptic  fianctions.  However,  Dr.  C.  Masai tis 
has  obse2"ved  that  if  is  parallel  to  a  principal  axis  of  C  then  w  Is 
expressible  in  terms  of  elliptic  functions. 
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n 

I  • 


Axi symmetric  Systems 

In  Section  6  we  assumed  that  all  characteristic  roots  of  the  matrix 
G  vspc  dJ-stinctj  &nd.  tiisi't  J*  V7s.s  psjrsUd  to  nono  of  ths  px’incipsl  sjccs 
of  C.  If  D  5^  E  =  F  and  +  i^)  0,  for  example,  then  (6.9)  takes 

the  degenerate  form 

d£_  =2T./r2  (7.1) 


We  obtain  a  similar  enuation  for  the  determination  of  X  if  we  assume 
that  D,  E,  F  are  distinct,  but  exactly  one  of  the  components  of  J* 
vanishes.  Now  for  small  j  r  |  the  analysis  starts  with  four  critical 
points,  one  of  which  must  be  a  saddle- point.  With  increasing  j  r  j 
we  pass  finally  to  two  critical  points,  just  as  in  the  more  general 
circumstances  considered  in  Section  6. 

If  D  =  E  =  F,  then  (6.6)  degenerates  to 


DJ*^/(X-D)^  =  2T*/  r  ^ 


(7.2) 


Wc  obtain  a  similar  equation,  regardless  of  the  nature  of  the 
characteristic  roots  of  C,  if  only  one  component  of  J*  is  non-zero, 
or  if  E  =  F  and  j  =  0.  A  complete  enumeration  of  the  possibilities 
has  no  especial  interest,  and  in  any  event,  all  types  that  can  arise 
have  been  mentioned  already.  Now  there  are  always  only  two  critical 
points . 


To  turn  to  the  most  important  of  these  degenerate  cases,  suppose 
our  liquid-solid  system  is  axi symmetric,  with  respect  to  the  x-axis. 
Much  more  can  now  be  said  about  form  of  the  velocity  potential.  Let 
us  introduce  cylindrical  polar  coordinates 


i© 

x=x,  y+iz=re 


If  we  let  n  be  the  radial  component  of  the  inward  normal  to  S 


L' 


then 


n  +  in 

V  *7 

J  “ 


=  n  e 


Since  now  (x  x  n)^  =  0,  we  observe  that 


cj  =  0 
X 

satisfies  the  relevant  parts  of  equations  (2-5)  (2.6). 

(x  X  n)  -  i(x  X  n)  =  (n  x  -  n  r)  e^^ 

—  z  —  —  y  r  X  ' 


(7.5) 


Since 


a  (x)  -  ia  (x)  =  '|f(x,r)  e 
z  -  y  - 


iw 


will  satisfy  (2.5)  and  (2.6)  if 


(7.4) 


1 

at 

t 

ax^  ar^ 

r 

ar 

■  2  ■ 
r 

at 

n  +  n 

ax  '  ‘ 

at 

n  X 

—  n  T* 

“x* 

ar 

**r“ 

in  S* 


nn  f!* 

““  ”L 


(7.5) 


(7 

\  t  •  / 


where  S*  is  a  cross  section  of  V.^  in  6iny  plane  0  =  constant,  and  C*  is 

Li  Li 

its  boundary.  Obviously 

t(x)  =  9/2it  (7.7) 

satisfies  (2.9)  and  (2.1o)  and  increases  by  unity  for  each  positively- 

directed  circuit  of  a  circle  x  =  x  ,  r  =  r  . 

o  o 

To  write  the  equations  of  motion  we  would  require 


rp  ^ 

2T*  =  w^Cj^  w  =  J 


^  (  V  w  •£)  d.V 


By  (7-3)  and  (7-4)  this  becomes 

2T*  =  jtp^  (w^^  +  w^)  r  To*  f  +  C  + 
. ■  ■ 


r  drdx 


By  means  of  (7-5)  and  (7.6),  and  Gauss'  theorem  this  can  also  be 
written  in  the  more  convenient  form  for  calculation 


(7.8) 


2T*  = 


np^  (w*^  +  w  )  /  rt  (n  X 
y  z'  J  ^  r 

L> 


n  r)  ds 

X 
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(7.9) 


where  s  is  arc-length  along  C*.  Clearly  the  matrix  C  is  proportional 

Li  Li 

to 

0  0  0 

/  \ 

(  0  .  0  ) 

'  0  0  1  ' 

We  also  need  J*.  Since  for  (7-7) 

-2rt  (x  X  Vt)  =  (-1,  X  cos  o/r,  x  sin  S/r) 


then 


J*  =  -  y  y  X  X  V  'T  dv  yields 


-  J*  =  (p^  jj  g*  r  drdx,  0,  0)  (7-10) 


It  should  also  be  remarked  that  the  transverse  moments  of  inertia  of 


the  solid  Eg  =  Fg,  i.e. 


Dg  0  0 

0  Eg  0 

0  0  E. 


\ 


Thus  C  =  C,,  +  C  will  also  be  diagonal,  and  the  last  two  diagonal 
Id  1-1 

elements  will  be  equal. 

Uow  (6.2)  and  (6.5)  become 

D  w^  +  E  (w  w  ^)  =  2T* 

X  y  z 


(7.11) 


^  2  _2  /  2  2 
(D  w  +  r  j)  +  E  (w  +  w  )  =  K 


f  r-r  ^  \ 

W  .  -L-i  } 

where  j  is  the  x-coraponent  of  J*  in  (7.8).  The  intersections  of  these 
ellipsoids  of  revolutions  are  circles 

w  =  constant 
X 

2  2  „  2w„  _  ^’2. 

w  +  w  =  -  u  w  ;/n  =  a 

y  z  '  X 


?5 


2 


2 


Thus  the  solutions  of  the  equations  of  motion  (6.1)  become 


V 

y 


constant 


Iv  (t-t  ) 
e  o 


vhere  the  rate  of  precession 

V  =  (  1  -  d/e)v^  -  r  j/E 


(7-13) 


(7.14) 


For  practical  purposes  one  would  certainly  be  primarily  interested 
in  axisymmetrlc  systems.  However,  it  should  be  remarked  that  small 
errors  will  occur  in  machining  models  intended  to  be  axisymmetrlc.  If 
many  models  are  to  be  constructed,  one  may  also  deliberately  abandon 
nearly  perfect  dynamic  axi symmetry  because  the  effort  to  achieve  it 
does  not  result  in  sufficient  Improvement  in  the  performance  of  the  model. 
Therefore  the  more  general  analysis  of  Section  6  may  also  have  some 
bearing  on  practical  applications.  For  systems  that  are  not  quite 
axisymmetrlc  one  would  expect  the  full  sequence  of  possibilities 
depicted  in  Figure  2  to  occur.  In  Figures  2a  and  2b,  however,  the 
areas  occupied  by  closed  integral  curves  surrounding  and  m^  would 
presxjmably  cover  most  of  the  surface  of  the  energy  ellipsoid,  and  the 
remaining  sets  of  Integral  curves  should  be  cramped  into  a  relatively 
small  fraction  of  the  surface  area. 
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8.  TOPOLOGICALLY  COMPLICATED  CAVITIES 


Let  us  motivate  the  following  discussion  by  beginning  with  a 
special  configuration.  Consider  a  solid  with  a  toroidal  cavity  V 

Ij 

bounded  by  the  concentric  cylinders  r  =  r^  +  r  and  the  parallel 
planes  x  =  +  x^.  To  inodify  the  motion  of  the  solid- liquid  system  let 
us  insert  partitions  on  the  planes  9  =  2nm/n,  1  »  m  =  n.  This  will 
merely  subdivide  the  cavity  into  n  topologically  spherical  regions^ 
which  creates  no  dynamical  novelty.  Suppose,  however,  that  all  of  the 
partitions  are  perforated.  If  each  partition  contains  exactly  one  hole, 
will  be  topologically  toroidal  (l.e.  homeoraorphlc  to  the  interior  of 
a  torus ) .  But  if  we  suppose  that  some  of  the  partitions  contain  more  than 
one  hole  the  topological  structure  of  will  become  more  complicated 
than  the  interior  of  a  torus.  Since  the  following  discussion  will  not 
be  concerned  with  the  exact  number  of  perforations,  we  shall  emphasize 
its  generality  by  suggesting  that  the  partitions  could  even  be 
imagined  to  be  made  of  finely  woven  wire  mesh,  to  impart  an  extremely 
complicated  topological  structure  to  Vj^. 


Now  let  us  indicate  how  the  treatment  of  the  toroidal  cavity  can  be 
adapted  to  the  case  of  a  very  general  cavity  not  necessarily 
constructed  by  the  process  described  in  the  preceding  paragraph.  In  our 
earlier  discussion,  the  topological  nature  of  V^.  asserted  itself 
through  applications  of  Gauss'  theorem  to  various  volume  integrals 
taken  over  Let  us  suppose  that  is  bounded  by  a  two-sided  surface 

of  genus  N.  In  other  words,  can  be  considered  to  be  topologically 
equivalent  to  the  interior  of  a  sphere  with  N  handles.  Then,  by  the 


insertion  or  N  partitions  S*  we  can  make  v'  into  a  topological  sphere 

J  L 

V' .  With  each  partition  S*  we  associate  a  velocity  potential  fianction 
L  J 

T.(x)  which  produces  unit  circulation  on  a  closed  path  in  V'  from  the 

tJ 

(arbitrarily  chosen)  initial  side  to  the  final  side  of  S*.  In  fact, 

J 

T.  Will  be  single-valued  in  the  cavity  of  genus  N-1  produced  by 

tJ 

Inserting  only  the  partition  S*  but  none  of  the  others.  Now  replace 
/  Q  0\  f  '^”1 


<I>  (x,t)  =  X  .  u(t)  +  w(t)  .0  (x)  +  T  r  T  (x)  (8.1) 

J  J 
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where  r .  is  a  constant  circulation  associated  with  t..  The  boundary 
J  J 

conditions  tor  a  (x)  are  as  before,  while,  analogously  to  (2.10) 


Sr^/dn  =  0 
J 


on  S, 


^6.2) 


The  kinetic  enerstv  of  the  solid-liquid  system  becomes 


2T  =  Mu  + 


T 

wCw  +  va,  ,  r  ^  r  . 


(B.5) 


where 


=  /  V  •  Vt  dV 

'  L 


(8.4) 


is  a  constant  positive-definite  matrix.  As  before,  the  linear 
momentum  is 


-  dT/du  =  M  u  (8.5) 

and  the  angular  momentum 

=  dT/dw  + 

where  the  axis  vector 

"  /v^  Pf.  -  ^ 

is  associated  with  t.. 

J 

For  motion  subject  to  no  external  moment  we  again  obtain  the 
energy  integral 


(8.6) 


T^ 

W  UW  = 

and  the  angular  momentum  integral 

(Cw  +  y 

The  search  for  the  critical  points, 
/A  Ai  +0 

w  Ji  (J  ^  y  t/W 


J  -J 


w  ,  of 
— c 


(8.8) 

=  K  (8.9) 

the  function  K  on  the  energy 


or 


w^.-(c-ur  iTjJ. 


(8..11) 


When  we  substitute  (8,ll)  into  (8,8)  we  obtain  as  an  analog  of  (6.6) 


(8.12) 


where 


P(X)  =  (C  -  XI)”^  C(C  -  XI)“^ 


(8.13) 


The  analysis  of  the  nature  of  the  integral  coi'-ves  as  a  function  of 
the  N  parameters  P  cein  be  carried  out  along  the  following  lines. 

Suppose  that  n(  s  N)  of  the  vectors  J*  are  linearly  independent,  where 

<  <  <  < 

1  =  n  =  3.  Let  K  ,  for  1  s  u  s  n,  be  an  orthonormal  basis  for  the  set 
-a’ 

J*.  Then  we  must  have 

""  J 


•i,  -  1.  Sa-  i!p  ■  °>  “  P 


(8.1L) 


Also,  there  must  exist  an  N  x  n  matrix  G  of  constant  elements,  and  of 

jOC 

rank  n,  such  that 


J*  =  y  ?  G  H 
-j  ^1  jo-a 


/  Q  nr-' 
j-p; 


Now 


N 


^  1  J  ^  1 


n  N 

y  /  y  ^  i~.  \  rj 

^  1  j  -o 


Let  W  be  any  n-diraensional  unit  vector,  i.e. 

a 


y  ”  =1 


(8.16) 


There  is,  of  course,  an  n-1  parameter  family  of  W^.  Then  for  any  P 
the  system  of  linear  equations 


n 


r  G,  =  r  w 

J  ja  a 


(6.17) 
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has  an  (N-n) -parameter  family  of  solutions  r  ^  .  The  general  solution 
of  (6.17)  is  of  the  form 


N-n 

f  •‘•I  K  y  if  r  0 


N-n 


■  j 


V 

4^ 


Je 


if 


=  0 


(8.18) 


where  y  . 

Jo 


is  a  particular  solution  of 


Jo 


and  if  N  >  n,  then  7.  are  N-n  linearly  independent  solutions  of  the 

J£ 

associated  homogeneous  equations,  and  A  are  N-n  arbitrary  constants. 

■  '  e 

Now 


N 

^  Z  ”  w  H 

1  J^j  1  a  -a 


and  (8.12)  takes  the  form 


n  T  n 

Z  W  H  P(X)  S  W  H  =  2T*/  r 


if  r  /  0 


(8.19) 


(8.20) 


which  is  now  more  closely  analogous  to  (6.6).  The  kinetic  energy  due  to 
circulation  is 


'•5  1 V,  '"i  >■  j  ■ 


r  ^  if  r  7^  0 

d^  if  r  =  0 


/  «  DT  ^ 


where 


N-n 


c  = 


N-n 


0.5yA  (7  +y  A  7^)(/,  ..... 

^  ij  ^ 'io  ^  €  ie'  Jo  y  ’  jy 


d"^  =  0 . 5  y  A  T  Ay  y 


N-n 


A_  y.^) 

■■  -jy' 

(8.22) 


€  re 


A^y 

y  jy 


Thus  T*/  r  G  would  be  the  ratio  of  kinetic  energy  due  to  w  to  that  due 

2  2 

to  circulation  when  r  ^  0.  Note  that  if  N  >  n,  then  c  and  d  may  vary 

with  the  choice  of  A  . 

c 
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If  we  make  particular  choices  of  and  c  ,  then  in  general  we  have 

exactly  the  relation  between  F  and  the  integral  curves  on  the  energy 

ellipsoid  that  is  described  in  Section  6.  If  N-n  =  2,  (8.22)  has  an 

(N-n-l)  -  parameter  family  of  solutions  .  This  simply  means  that  the 

same  set  of  motions  of  the  system  can  be  realized  with  an  (N-n-l)  - 

parameter  set  of  choices  of  the  circulations  F,. 

J 

Since  equation  (8.20)  is  the  crucial  element  in  the  discussion  of 
possible  motions  of  the  solid- liquid  system,  then  the  categorization  of 
motions  should  clearly  be  based  on  whether  F  =  0  or  F  /  0,  and  then,  in 
the  latter  case,  on  the  number  of  parameters,  n-1,  required  to  determine 
W^.  Thus  there  will  be  four  major  types  of  motion.  The  distinctions 
between  them  could  conceivably  be  visualized  and  clarified  by  describing 
some  of  their  properties,  such  as  (l)  the  locus  of  the  centers 

w  =  -  C"^  y  F  J*  =  -  F  y  W  C"^  H  (8.25) 

of  the  angular  momentum  ellipsoids  (8.9)  as  a  function  of  F  and  W^; 

(2)  the  possible  steady  states  of  rotation  for  fixed  F  and  (if  possible) 
variable  W^J  and  (3)  the  limiting  steady  states  for  |  F  j  =00. 

CASE  I.  If  F  =  0,  then  by  (8.23)  the  center  of  the  angular  momentum 

ellipsoid  is  at  the  origin.  If  N  >  n  there  may  actually  be  circuits  with 

circulation  F  0,  in  accordance  with  the  second  part  of  (8.I8).  By  (8.I9) 
J 

for  F  =0  and  (8.IO)  \  must  be  a  characteristic  root  of  C,  and  w^  a 
characteristic  vector.  This  leads  to  the  familiar  rigid  body  and  spherical 
cavity  type  of  motion. 

CASE  II.  If  F  ^0  and  n  =  1,  the  results  of  Section  6  are  applicable  word 
for  word.  By  (8.25)  the  centers  of  the  momentiom  ellipsoids  are  on  a 
straight  line  through  the  origin  pacrallel  to  C  For  fixed  F  there 

are  from  two  to  six  possible  steady  states  of  rotation,  depending  on  the 
magnitude  of  F  .  For  |  •”  |  =  oo  there  are  two  steady  states  of  rotation 
at  the  ends  of  the  diameter  of  the  energy  ellipsoid  parallel  to  the  vector 
C  .  Such  motions  occur,  in  particular,  for  toroidal  cavities. 
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CASE  III.  If  r“  ^  0  and  n  =  2,  by  (8.23)  the  centers  of  the  angular 
momentum  ellipsoids  are  on  a  plane  through  the  origin.  Let  us  set 
=  cos  0,  =  sin  ©.  For  fixed  T  and  0  there  are  from  two  to  six 

steady  states  of  rotation,  depending  on  the  value  of  r  .  If  ©  varies 
while  r  remains  fixed  the  critical  points 

w^  =  -  r  [  c  -  X(©,  r  )  I j  (cos  9  +  sin  ©  H2) 

will  trace  a  set  of  curves  on  the  energy  ellipsoid.  In  accordance 
with  the  results  obtained  at  the  end  of  Section  6,  for  |  T  j  =  co  and 
fixed  ©  the  two  corresponding  steady  states  of  rotation  will  be  at  the 
ends  of  the  diameter  parallel  to  V W  H  =  cos  ©  H  +  sin  ©  H„.  In  other 
words,  for  |  I”  |  =00  the  possible  steady  states  are  on  the  intersection 

of  the  energy  ellipsoid  and  the  plane  ^  w  =  0. 

CASE  IV.  If  r  p  and  n  =  5j  then  the  centers  of  the  angular  momentum 
ellipsoids  can  be  anywhere  in  w- space.  Let  us  set  =  cos  ^  cos  ©, 

Wg  =  cos  ^  sin  ©  ,  Wj  =  sin  ©.  For  fixed  T  ,  ^  and  ©,  there  will  be 
from  two  to  six  steady  states  of  rotation,  depending  on  the  value  of  r  . 
If  ^  and  ©  vary  Independently  while  P  remains  fixed,  the  critical 
points 


will  trace  out  regions  on  the  energy  ellipsoid.  In  accordance  with  the 
limits  calculated  at  the  end  of  Section  6,  for  fixed  ^  and  ©  the  two  steady 
states  of  rotation  for  |  T  j  =00  will  be  at  the  ends  of  the  diameter 
of  the  energy  ellipsoid  parallel  to 

cos  ^  (cos  ©  +  sin  9  H2^  ^  ^5 

If  we  let  ^  and  9  range  over  all  permissible  values,  we  obtain  the  entire 
surface  of  the  energy  ellipsoid. 

Let  us  conclude  by  reiterating  that  Case  I  Includes  the  topologically 
spherical  cavity  (genus  zero),  eind  Case  II  the  toroidal  cavity  (genus  one). 
It  certainly  seems  reasonable  to  conjecture  that  Cases  III  and  IV,  respec¬ 
tively,  correspond  to  cavities,  or  at  least  to  some  cavities,  of  genera 


two  and  three,  respectively.  In  a  cavity  of  genus  N  it  is  possible  to 

assign  arbitrarily  N  independent  circulations  F  .  For  N  ^  4  our 

J 

conjecture  would  imply  that  increases  in  genus  do  not  lead  to  new  types 
of  dynamic  behavior,  but  merely  present  a  greater  variety  of  choices  of 
parameters  (  F  )  to  simulate  the  behavior  of  bodies  with  liquid  filled 

tJ 

cavities  of  genera  less  than  four.  If  there  were  only  some  mechanism 
for  randomly  exciting  and  varying  strong  circulations  in  a  body  with  a 
complicated,  liquid  filled  labyrinth,  one  might  speculate  that  then  by 
virtue  of  the  possibility  of  suddenly  inducing  degeneracies  of  the  sort 
discussed  in  Section  7>  thereby  switching  from  motion  of  one  type  to 
another,  the  behavior  of  the  liquid-solid  system  could  become  highly 
erratic  and  unstable. 


J.  H.  GIESE 
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